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Abstract. We study the magneto-rotational instability of an incompressible flow which rotates with angular 
velocity Q(r) = a + b/r 2 where r is the radius and a an d b are constants. We fin d that an applied magnetic 
field destabilises the flow, in agreement with the results of Riidiger & Zhang (2001). We extend the investigation 
in the region of parameter space which is Rayleigh stable. We also study the instability at values of magnetic 
Prandtl number which are much larger and smaller than Riidiger & Zhang. Large magnetic Prandtl numbers are 
motivated by their possible relevance in the central region of galaxies (Kulsrud & Anderson 1992) . In this regime 
we find that increasing the magnetic Prandtl number greatly enhances the instability; the stability boundary drops 
below the Rayleigh line and tends toward the solid body rotation line. Very small magnetic Prandtl numbers are 
motivated by the current MHD dynamo experiments performed using liquid sodium and gallium. Our finding 
in this regime confirms Riidiger & Zhang's conjecture that the linear magneto-rotational instability and the 
nonlinear hydrodynamical instability (Richard & Zahn 1999) take place at Reynolds numbers of the same order 
of magnitude. 
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1. Motivation 



It is thought that turbulence in accretion discs arises 
from a magneto-rotational instability (MRI), where a 
magnetic field destabilises a rotating velocity field which 
decreases outwa rdly. This instability was d iscovered by 
|Yolikhov (1959) and Chandrasckhar (1961) when study- 
ing the motion of a fluid between two concentric cylin- 
ders (Taylor-Couette flow). It was only years later, when 
it was realized that Couette flow can be interpreted as 
model of Keplerian flow, that the implications for astro- 
physics were fully appreciated (Balbus & Hawley 1991). 
There have been many studies of the MRI recently, from 
numerical simulations of accretion discs (Brandenburg et 
al. 1995) to nonlinear calculations in spherical geometry 
(Drecker et al. 2000). In particular, in a recent paper 
Riidiger & Zhang (2001) analysed the linear stability of 
hydromagnetic Couette flow and showed that, if the mag- 
netic Prandtl number is less than unity, azimuthal Couette 
flow is more easily destabilised with a magnetic field than 
without. They found that the instability extends into the 
region of parameter space which is Rayleigh-stable with- 
out a magnetic field, which is important since Keplerian 
rotation is Rayleigh-stable. Riidiger & Zhang were also 
able to study the MRI instability at magnetic Prandtl 
numbers Pm as small as 0.001, towards the limit relevant 
to liquid sodium and gallium {Pm ss 10~ 5 ) which are used 
in current MHD dynamo experiments (Tilgner 2000). 



Th e aim of this paper to is to extend the investiga- 

We explore the instabil- 



tion of IRiidigcr fc Zhang (2001 



ity as a function of the speed of the outer cylinder in the 
Rayleigh-stable region, which is the parameter space of 
astrophysical interest (Riidiger & Zhang considered only 
one nonzero ratio of outer to inner cylinder's rotation). We 
also determine the effect of changing the magnetic Prandtl 
number, extending the range studied by Riidiger & Zhang. 



2. Model 

We consider an incompressible fluid contained between 
two coaxial cylinders of inner radius R\ and outer radius 
i?2 which rotate at prescribed angular velocities Sli and 
0,2- We use cylindrical coordinates (r,<fi,z) and assume 
that a uniform magnetic field of strength Bq is applied in 
the axial z direction. At small angular velocities the flow 
is purely azimuthal (circular-Couette flow) and has mag- 
nitude Vo(r) = ar + b/r which corresponds to the rotation 
law, 



(1) 



The constants a and b are determined by the no-slip 
boundary conditions at the cylinder's walls, so we have, 



a = Qi 



{R\/Rj - n 2 M) 

{l-Rl/Rl) ■ 



b = Q 1 Rf 



2 (1 



n 2 M) 



{i-Ri/Rir 
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At sufficiently high angular velocities Couette flow be- 
comes unstable. The resulting velocity and magnetic fields 
V(r, (f>, z, t) and B(r, <f), z, t) are determined by the MHD 
equations which we write in dimensionless form as, 



dV 
~dt 

dB 

~dt 



(V-V)V = -Vp 



V 2 V 



Q_ 

Pm 



Pm 



V 2 B + V A (V A B), 



(VAB)A B, (3) 



(4) 



V • B = 0, V • V = 0, 



(5) 



where p is the pressure. In writing (^|) and (^) we used 
5 = R 2 — i?i as unit of length, S 2 jv as unit of time and 
Bq as unit of magnetic field, where v is the kinematic 
viscosity. The governing dimensionless parameters of the 
problem are the Reynolds numbers Re\ , Re 2 , radius ratio 
rj, and rotation ratio p, 



Rei = 



R^S 



Re 2 



R 2 n 2 6 



V 



R2 



/' : 



^2 



(6) 



together with Chandrasekhar's number Q and the mag- 
netic Prandtl number Pm. 



Q 



pv X 



(7) 



where p is the density, a is the electrical conductivity, 
A = l/(apo) is the magnetic diffusivity and po is the per- 
meability. 

We solve the MHD equations by direct time stepping 
from a small seeding initial condition and determine the 
stability of the Couette solution in different regions of pa- 
rameter space. Our numerical method for 3D nonlinear 



flow is detailed in Willis & Barenghi (2002). For this lin- 
ear study it suffices to say that the formulation is based 
on representing V and B with suitable potentials which 
are spectrally expanded over Fourier modes exp[iTO(/>+ia:z] 
in the azimuthal and axial directions and over Chebyshev 
polynomials in the radial direction. The axial wavelength 
of the disturbance is 2-k /a. We assume no-slip boundary 
conditions for V, and electrically insulating boundaries 
conditions for Z?, 



dB r 
dz 



dr 



1 dB z 



dB 4 
dz 



(8) 



where B m (r) denotes the modified Bessel functions, 
I m (ar) on R\, and K m (ar) on R 2 . The time stepping is 
based on a combination of second order accurate Crank- 
Nicolson and Adams-Bashforth methods. The resulting 
computer code was tested against published results with 
and without magnetic field (Chandrasekhar 1961; Roberts 
1964; Marcus 1984; Jones 1985; Barenghi 1991). 

3. Results 

It is well known that, in the absence of magnetic field, 
inviscid Couette flow is linearly stable provided that the 
celebrated Rayleigh criterion (p > rj 2 ) is satisfied. 




Fig. 1. rj = 0.5, Pm = 1. Critical Reynolds number of the 
inner cylinder (i?ei c ) versus applied magnetic field (^/Q); 
two cases are presented: outer cylinder fixed (p — 0) and 
on the Rayleigh line (p = rj 2 ). 



Figure 1 shows the result of our calculations for the 
stability of dissipative Couette flow to axisymmetric per- 
turbations, for radius ratio rj = 0.5 and magnetic Prandtl 
number Pm = 1, as a function of the applied magnetic 
field. We plot the result in terms of y/Q rather than 
Q in order to make direct comparison with the work of 
Riidigcr fc Zhang (2001) . The first curve refers to the case 
in which the outer cylinder is fixed, p = 0, which is the 
most studied case in the fluid dynamics literature. It is 
apparent that the presence of a magnetic field makes the 
flow more unstable. The critical Reynolds number, which 
is Rei c — 68.2 for Q = 0, decreases with increasing Q and 
has a minimum at Q = 39. The most unstable mode is 
m = over the range for Q in Fig. 1. From here we con- 
sider axisymmetric disturbances only. The critical axial 
wavenumber a decreases significantly from 3.1 to 1.7 over 
the range, and varies like 1/VQ thereafter. This stiffening 
eventually restabilises the flow and Re\ c increases like \/Q 
for strong fields. 

The initial destabilisation is consistent with the find- 



ing of Riidiger fc Zhang (2001) ; the small difference be- 
tween their Re\ c and ours is certainly due to the different 
boundary conditions for B (they assumed pseudo-vacuum 
conditions and we assume insulating conditions). The sec- 
ond curve of Fig. 1 refers to the case p — rj 2 (the Rayleigh 
line), which separates stable and unstable regions in the 
absence of a magnetic field. The curve well illustrates the 
striking destabilising effect of the magnetic field. 

Figure 2 shows stability boundaries in the Re\ vs Re 2 
plane for different values of the imposed magnetic field. 
The axial wavenumber a does not vary a great deal along 
the boundaries, its dependence being principally deter- 
mined by the strength of the field Q. 

It is apparent that even a small value of Q is enough 
to make the boundary cross the Rayleigh line (the upper 
dotted curve in the figure). The destabilising effect of the 
magnetic field is so large that the stability boundary drops 
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Fig. 2. i] = 0.5, Pm = 1. Critical Reynolds number of the 
inner cylinder (Rei c ) versus Reynolds number of the outer 
cylinder (Pea) at different values of applied magnetic field 
(Q). The upper dotted line is the Rayleigh criterion. The 
lower dotted line is solid-body rotation. 




Fig. 3. i] = 0.5, Q = 10. Critical Reynolds number of the 
inner cylinder (Pe lc ) versus Reynolds number of the outer 
cylinder (Re2) at different values of magnetic Prandtl 
number (Pm). The upper dotted line is the Rayleigh cri- 
terion. The lower dotted line is solid-body rotation. 



toward the region of solid body rotation (Sli = Q2 or 
[i = 1), which is the lower dotted line. However, if the 
applied field is strong enough the flow can be restabilised, 
in accordance with the results of Fig. 1. 

The destabilisation becomes more important the larger 
the magnetic Prandtl number Pm. Figure 3 shows results 
for a moderate value of applied magnetic field, Q — 10, at 
increasing values of Pm. It is apparent that the stability 
boundary drops much below the Rayleigh line [i = rj 2 
and, for large enough Pm, becomes asymptotic to the solid 
body rotation line fi = 1. 

Finally, in Fig. 4 we find that on the Rayleigh line the 
stability follows a power law, 




Fig. 4. T) = 0.5. Critical Reynolds number of the inner 
cylinder (Re\ c ) versus magnetic Prandtl number (Pm) on 
the Rayleigh line, /i = rj 2 . 



over a surprising range of Pm. The error in j3 is based on a 
fit through all points calculated. For 7/ = 0.5 the minimum 
Rei c occurs at Q = 86 and a = 2.1. To verify this result 



we use the WKB analysis of Ji et al. (2001). The local 
dispersion relation they find for the Taylor-Couette flow 
is identical to the one derived for accretion discs in the 
incompressible limit (Sano & Miyama 1999). From the 
dispersion relation the necessary and sufficient condition 
for stability to axisymmetric perturbations is, 

(Pm + S 2 ) 2 (l + e 2 ) + 2CPm 2 - 2(2 - QRm 2 S 2 > 0, (10) 

where S is the ratio of diffusive to Alfven timescales, Rm is 
the magnetic Reynolds number, and ( = l/(rf2) d r (r 2 il). 
On the Rayleigh line, from (jl|) and f2(r) takes the 
simple form, £1 = b/r 2 , and so £ = at all r. We said ear- 
lier that in Fig. 2 the wavenumber changes little along the 
boundary, as its dependence is principally on the strength 
of the imposed field. The wavenumber parameter e of 
Ji ct al. (2001) then depends only on Q, so e = e(Q). 
Ignoring factors that depend only on e in the substitu- 
tions, 



S = PmQ, Rm—PmRe, 



( |10| ) leads to the condition for stability, 

< m, 



(11) 



(12) 



which confirms our result ([)]). This relation ( |l2| ) should 
hold for all Pm as we have not been required to take any 
limits. The function A(Q) has a minimum for some Q 
and a(Q). For rj = 0.5 we found the minimum occurs at 
Q = 86, a = 2.1. 



Re lc ocPm~' j , (3 = -0.500 ± 0.002, 



(9) 



Whilst the linear result (12) should hold for very small 
Pm, from the work of Wendt (1933) and Taylor (1936) 
a turbulent (nonlinear) instability is observed at large 
Reynolds numbers (Pei > 2.5 x 10 5 ) with this radius ra- 
tio. Extrapolating down one magnitude to Pm = 10~ 5 for 
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laboratory fluids gives Re\ c w 2 x 10 5 , a similar value, for 
the magneto-rotational (linear) instability. 

A radial truncation of 12 Chebyshev modes was found 
sufficient for all calculations. In particular, for the calcu- 
lation at Pm = 1CP 4 in Fig. 4, convergence was tested by 
increasing the truncation and also decreasing the timestep. 
The fractional numerical error in Re\ c is estimated at ap- 
proximately 1CP 6 . 
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4. Discussion 

Our calculations show that many rotation laws of the form 
n(r) = a + b/r 2 which are hydrodynamically stable (that 
is to say, they satisfy the Rayleigh criterion) become lin- 
early unstable when a magnetic field is applied. Our re- 



sults confirm the finding of Riidiger & Zhang (2001) and 



extend them in the Rayleigh stable region. 

We have determined the instability at magnetic 
Prandtl numbers Pm one order of magnitude smaller than 
Riidiger & Zhang's, towards the small magnetic Prandtl 
number limit, which is relevant to possible MHD dynamo 
experiments with liquid sodium and gallium. Although the 
power law Re\ c oc Pm~ 5 that we find on the Rayleigh line 
(fi = T] 2 ) is slightly different from theirs (Re\ c oc PmT ' 65 
on /i = 4), it confirms their conjecture that the nonlinear 



instability found by Richards fc Zahn (1999) and the MRI 



are likely to occur at Reynolds numbers of the same order 
of magnitude. 

We also find that the flow becomes particularly unsta- 
ble if the magnetic Prandtl number is greater than unity. 
The instability boundary in the Re\ vs Re2 plane rapidly 
tends towards the solid body rotation line. This enhanced 
instability for large Pm is consistent with earlier results of 



Kurzweg (1963) . His boundary conditions were selected 



such as to avoid mathematical difficulties but for small 



Pm agreed well with the results of Chandrasckhar (1961) 
The significance of the instability in this case is linked to 
the possibility (Kulsrud & Anderson 1992; Brandenburg 
2001) that large values of Pm exist in central regions of 
galaxies. 
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